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Abstract 

We present simulations of congested traffic in circular and open systems with a non- 
local, gas-kinetic-based traffic model and a novel car-following model. The model parame- 
ters are all intuitive and can be easily calibrated. Micro- and macrosimulations with these 
models for identical vehicles on a single lane produce the same traffic states, which also 
qualitatively agree with empirical traffic observations. Moreover, the phase diagrams of 
traffic states in the presence of bottlenecks for the microscopic car-following model and the 
macroscopic gas-kinetic-based model almost agree. In both cases, we found metastable 
regimes, spatially coexistent states, and a small region of tristability. The distinction of 
different types of vehicles (cars and long vehicles) yields additional insight and allows 
to reproduce empirical data even more realistically, including the observed fluctuation 
properties of traffic flows like the wide scattering of congested traffic data. 

Finally, as an alternative to the gas-kinetic approach, we propose a new scheme for 
deriving non-local macroscopic traffic models from given microscopic car-following models. 
Assuming identical (macroscopic) initial and boundary conditions, we show that there are 
microscopic models for which the corresponding macroscopic version displays an almost 
identical dynamics. This enables us to combine micro- and macrosimulations of road 
sections by simple algorithms, and even to simulate them simultaneously. 

Keywords: Congested traffic flow, vehicle queues, phase diagram, micro-macro link, 
car-following model 

1 Introduction 

In recent years, the community of scientists engaged in traffic modelling is rapidly growing (for 
overviews see |2], 0]). This is not only due to practical implications for optimizing freeway 
traffic ||, but also because of the observed non-equilibrium phase transitions |5], H and the 
various non-linear dynamical phenomena like the formation of traffic jams J?J, stop-and-go 
traffic |], H 10], and "synchronized" congested traffic |JTT], [5], |12]. It seems that all forms of 



congested traffic have almost universal properties which are largely independent of the initial 
conditions and the spatially averaged density, like the characteristic outflow Q out from traffic 
jams of about 1800 ± 300 vehicles per hour and lane or their typical dissolution velocity C of 



about — 15 ± 5 kilometers per hour |13|]. This universality arises from the highly correlated 
state of motion produced by traffic congestions. In particular, the outflow Q ou t is related to the 
time gap between successive departures from the traffic jam |T4|, |T3|. Therefore, the outflow 



is practically independent of the initial conditions and the kind of congested traffic. As a 



consequence of the constant outflow, the propagation velocity C of jam fronts, given by the 
dissolution speed of traffic jams, is nearly constant as well. 

From a physics point of view, a model for a real system (here: for unidirectional freeway 
traffic) should be as simple as possible, but not simpler (i.e., parameter dependencies that are 
known to exist, should be reflected by the model). The parameters of the model should be 
intuitive, easy to calibrate, and the corresponding values should be realistic. Moreover, the 
traffic model should reproduce all observed localized and extended traffic states [TBj , including 
synchronized traffic and the wide scattering of congested traffic data [D]]. Furthermore, the 
observed hysteresis effects [[H], flSfl , complex states [11], |l9| , and the existence of self-organized 



quantities like the constant propagation velocity of stop-and-go waves or the outflow from a 
traffic jam QXO | should be reproduced. Finally, the dynamics must not lead to vehicle collisions or 
exceed the maximum vehicle density, and the model should allow for a fast numerical simulation. 

In this paper, we discuss some of these criteria for the recently proposed non-local, gas- 
kinetic-based traffic model (GKT model) ]T2|, ED[ and the microscopic intelligent-driver model 



(IDM) ||21|| . In particular, we discuss several different kinds of bottlenecks and present the 
phase diagram of congested traffic states for open, inhomogeneous systems p2f. We show that, 
for identical vehicles, the macroscopic dynamics is qualitatively the same for both models. Fi- 
nally, we investigate heterogeneous (mixed multi-class traffic composed of cars and lorries/long 
vehicles) and propose a natural explanation of the observed wide scattering of congested traffic 
data. 



2 Macroscopic Traffic Models 

In contrast to microscopic traffic models, which delineate the positions x a (t) and velocities 
v a (t) of all interacting vehicles a, macroscopic traffic models restrict to the description of the 
collective vehicle dynamics in terms of the spatial vehicle density p(x, t) and the average velocity 
V(x, t) as a function of the freeway location x and time t. One advantage of macroscopic traffic 
models is that they allow to simulate the traffic dynamics in several lanes by effective one- 
lane models considering a certain probability of overtaking. In the following, we will use this 
approach, although it is also possible to develop macroscopic traffic models for several lanes 
with explicit treatment of lane-changing, as mentioned later on. 

Being a consequence of the conservation of the number of vehicles, virtually all macroscopic 
traffic models are based on the continuity equation 

for the vehicle density p(x,t). The source term u(x,t)dx is the rate of vehicles entering or 
leaving the freeway at on- or off-ramp sections of length dx. The differences between the 
various existing macroscopic traffic models mainly concern the equation for the average vehicle 
velocity V(x,t). The Lighthill-Whitham model [H, H and its variants H, ||, |7|, H, H 



assume the equilibrium relation V(x,t) = V e (p(x,t)). For the description of emergent traffic 
jams and stop-and-go traffic, however, one needs a dynamic velocity equation. Most proposals 
can (in their continuous form) be summarized by the velocity equation 

dV y dV ldP l( V _ V ) (2 ) 

dt s_^S P ® x t 

Transport Term Pressure Term Relaxation Term 



Their main difference is the specification of the traffic pressure P, the relaxation time r, and 
the dynamic equilibrium velocity V e , which depends on the local vehicle density p. Notice 



that the Lighthill-Whitham model results in the limit r — > 0. Payne's [|30| and Papageorgiou's 
[3lH model is obtained for P(p) = [V — V c (p)]/(2t), with the "free" or "desired" average 
velocity V = V e (0). For dP/dp = — p/[2r(p + K)]dV e /dp, one ends up with Cremer's |3"2| 
model. In the model of Phillips ||33|| , there is P = p6, where 9 denotes the velocity variance. 
The model of Kuhne IB, PI, Kerner and Konhauser ITU, and Lee et al. |33j (KKKL model) 



results for P = p6o — r/dV/dx, where 6q is a positive constant and ry some viscosity coefficient. 
In comparison with a similar model by Whitham P5| , the contribution —rjdV/dx implies an 
additional viscosity term (r] / p) d 2 V / dx 2 . This is essential for smoothing shock fronts, which 
is desireable from empirical and numerical points of view (but causes some inconsistencies for 
large density or velocity gradients [[36], which can be avoided by a non- local macroscopic traffic 
model |7|). 

Many of the above mentioned models have proved their value in various applications, but 
they have been seriously criticized by Daganzo J36| because of theoretical inconsistencies. Apart 
from that, it is hard to decide which model is the best one. Whereas classical approaches 
focussed on reproducing the empirically observed velocity-density relation and the regime of 
unstable traffic flow, recent publications pointed out that it is more important to have traffic 
models which are able to describe the observed spectrum of non-linear phenomena and their 
characteristic properties [|7|, |K| [13], |2(| ^ . We think that it would be desireable to develop 
models that reproduce both aspects of empirical data. In the following, we will propose such a 



macroscopic traffic model, which is theoretically consistent as well J40 , [41 



2.1 The Non-Local, Gas-Kinetic-Based Traffic Model (GKT model) 

We derived macroscopic traffic equations from a gas-kinetic (Boltzmann-like) traffic model that 
was obtained from a simple microscopic car-following model under the approximation of quasi- 
instantaneous braking interactions [Q, |2(| (which will be circumvented by the micro-macro link 
suggested in Sec. §). Boltzmann-like traffic models for the spatio-temporal evolution of the 
so-called phase space density (= vehicle density p(x,t) times velocity distribution d(v;x,t)) 

45| , |46| , f|, [47|, £|8[ il, [5(J go back to Prigogine and coworkers (for an 
P2|). Paveri- Font ana improved Prigogine's model and cured some of its 
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Our gas-kinetic model can be viewed as an extension of Paveri-Fontana's 
equations to Enskog-like equations which are also valid in the regime of moderate and dense 
traffic, since they take into account the finite space requirements of vehicles. In particular, the 
equations become non-local to account for the fact that drivers react to the traffic situation 
ahead of them. This makes the derivation of the macroscopic traffic equations much more 
difficult than for the previous gas-kinetic models, the validity of which was limited to low vehicle 
densities. Therefore, our first calculations assumed "vehicular chaos" (i.e. uncorrelated vehicle 
velocities) and a gradient expansion in the vehicle density p(x,t), average velocity V(x,t) = 
J dv vd(v; x,t), and velocity variance 9(x,t) = J dv [v — V(x, t)] 2 d{y ; x,t) f4"If . In this way, we 
could derive necessary high-density corrections of the traffic pressure P(x,t) = p(x,t)9(x,t) 
(avoiding that vehicles would otherwise accelerate into jammed regions), and we could derive 
a plausible viscosity term, which had to be introduced in a phenomenological way in previous 
macroscopic traffic models. Meanwhile, we managed to carry out the calculations without a 

I, EH] , and even velocity correlations r among 



gradient expansion or other approximations |5T 



successive vehicles could be taken into account f53| (see Sec. |2/7]). We have also generalized our 
model to a multi-class multi-lane model 0, [53], which allows to account for different kinds of 
driver-vehicle units and lane changes in an explicit way. However, the results obtained up to 
now indicate that an efficient one-lane model (which averages over the neighboring lanes and 
the different driver- vehicle classes) can already capture the macroscopic vehicle dynamics in a 



semi-quantitative way [[53|, |4], Neverthess, in order to reproduce fluctuation effects like the 
wide scattering of congested flow-density data |], (□]], it is necessary to take the diversity of 
driver- vehicle types into account (see Sees. [2"T5| and |5TS|) . 

In the following, we will discuss the efficient one-lane version of our non-local, gas-kinetic- 
based traffic model for the case of uncorrelated vehicle velocities, but we will come back to 
correlation effects later on. The results of our excessive calculations can be represented in terms 
of the continuity equation (J]) and the velocity equation (|2|), but with a non-local, dynamical 
equilibrium velocity V e that depends on the vehicle density and other macroscopic quantities 
at the actual vehicle location x and at the advanced "interaction point" x'\ 

V e = V - r[l - p(p')]x(p')p'B(AV, S) . (3) 

v . ' 

Braking Term 

Apart from other things, this implies that V e implicitly depends on the (usually finite) variance 
9 and on the gradients of macroscopic variables like the density. A prime indicates that the 
corresponding variable is taken at the interaction point x' = (x + s) rather than at the actual 
position x. For simplicity, we have assumed the safe distance s to increase linearly with the 
average vehicle velocity V, 

S = 1 (—+TV) J (4) 

\ Pmax / 

where p max denotes the maximum vehicle density, T ~ 1.8 s the safe time headway, and 7 ~ 1 
an anticipation factor. 

The "Boltzmann factor" 

B(AV, S) = S { AVN{A V) + [1 + (AV) 2 ]E{ AV) } , (5) 

in which 

e -y 2 /2 

N(y) = (6) 
y An 

represents the normal distribution and 
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£'(,/) / dz Q -j= (7) 



the Gaussian error function, describes the dependence of the braking interaction on the effective 
dimensionless difference AV^ between the velocities at the actual position and the interaction 
point. The effective velocity difference is defined by 

V -V 

where 

S = 9 + 9'. (9) 

Hence, apart from the velocity difference to the interaction point, AV takes also into account 
the velocity variances 9 and 9' at the actual vehicle position x and the interaction point x', 
respectively. In spatially homogeneous traffic, we have B(0,S) = 5/2. In the limiting case 
AV 3> 0, where the preceding cars are much slower, we have a much stronger interaction 
B(AV, S) = (AV) 2 S. If, in contrast, the preceding cars are much faster (i.e. AV 0), we 
have B(AV,S) ~ 0. That is, since the distance is increasing, then, the vehicle at position x 
will not brake, even if its headway is smaller than the safe distance. 



Notice that, for 7=1, the macroscopic interaction term can be easily understood by 
the underlying microscopic dynamics of the gas-kinetic-based traffic model. If a vehicle at 
location x with velocity v is faster than one at x' with velocity w (i.e. Av = v — w > 0), 
it approaches the car in front within the time At = Axf iee /Av, where Aa:f ree = l/(p'x') is 
the average interaction-free distance per car. Then, if it cannot overtake immediately, which 
would happen with probability (1 —p'), it abruptly reduces the velocity by Av. The resulting 
ensemble-averaged deceleration is 

00 

(Av/ At) = -(l-p')xV/ d(Av) (Av) 2 d A (Av) . (10) 



If v and w are Gaussian distributed with averages V, V and variances 9, 9', respectively, 
the distribution function d A (Av) of the velocity difference Av is Gaussian distributed, with 
expectation value (V— V) = \/~S AV and variance S. Evaluating integral fllPl) yields (Av /At) = 
— (1 — p')x'p'B(AV, S), i.e., the macroscopic braking term in equation (||). 

Next, we will specify the pressure relation. Similar to gases, the calculations give 

P = p[9 + W(p)] (11) 



jq| . Again, 



= m = £ f 0i (12) 

1^1 1 P 

is the (weighted) average of the lane-specific velocity variances 9i, where the weights are the 
relative densities pi/ilp) in the I different lanes i. The monotonically decreasing, density- 
dependent function 

w = m - vf) = f bl (ytf - ff a v\ (13) 



can be used to correct for the variation of the average velocities V{ in the different lanes [see 
Figure |9|(a)], but it is often neglected. 

The pressure term describes the kinematic dispersion of the macroscopic velocity in inhomo- 
geneous traffic as a consequence of the finite velocity variance. For example, the macroscopic 
velocity in front of a small vehicle cluster will increase even if no individual vehicle accelerates, 
because the faster cars will leave the cluster behind. The kinematic dispersion also leads to a 
smooth density profile at the dissolution front between congested traffic and an empty road, as 
it occurs when a road blockage is removed (see Figure [I] and Ref. f37j). 

Insert Figure [7] about here. 

Notice that the total variance 

Q = (0 + W) (14) 

decreases with density and vanishes at the maximum vehicle density p max - Therefore, the 
density-gradient dPj dp of the traffic pressure P may decrease with growing density and increase 
with decreasing density. Hence, according to the pressure term (1/ p)dP/dx, it seems that 
vehicles would accelerate into denser regions and decelerate when driving into regions of lighter 
traffic. This tendency, however, is overcompensated for by the non-local interaction term 
Because of its non- locality, the braking term in relation (|3[) strongly increases with the 
density gradient. This can be more explicitly seen by applying the gradient expansions p' ~ 
p + s dp/dx + . . . etc. 



Next, in order to close the system of equations, one needs to specify the dynamics of the 
velocity variance 9. It turned out that, for a description of the presently known properties 
of traffic flows, we do not need to describe the variance by the associated partial differential 
equation. It seems sufficient to describe the variance by the equilibrium relation 

9 = Dt (15) 

of the dynamical variance equation |54 ] , where D denotes a velocity diffusion coefficient which 



is a function of the vehicle density p and the average velocity V. Empirical data suggest that 
the variance is a density-dependent fraction A(p) of the squared average velocity p0| , 



A(p) = A + AA 



tanh 



(17) 



6 = A(p)V 2 . (16) 

This guarantees that the velocity variance will vanish if the average velocity goes to zero, as 
required, but it will be finite otherwise. It turns out that the variance prefactor A is higher 
in congested traffic than in free traffic [Fig. 0(a)]. The empirical data can be approximated by 
the Fermi function 

' P-Pc 

where A and A + 2AA are about the variance prefactors for free and congested traffic, respec- 
tively, p c is roughly of the order of the critical density for the transition from free to congested 
traffic, and Ap denotes the width of the transition region. 

Insert Figure ^ about here. 

Finally, we will specify the "effective cross section" (1 — p)x = (1 —p)/p, assuming that 
the average time headway of vehicles corresponds to the safe time headway T, if the vehicle 
density is high (i.e. the freeway space is almost used up completely by the vehicular space 
requirements). This specification means that, in dense and homogeneous traffic situations with 
V = V and 6' = 9, we should have s = 1/p = [l/p m ax + 7V e (p)], where we insert the equilibrium 
solution fllDp. This gives us the relation 

W - v{p)\x{p) 



^4(Pmax)(l ~ P/P max ) 

which makes also sense in the low-density limit p — > 0, since it implies \ — > 1 and p — > 1. 



2.2 Properties of the Non-Local, Gas-Kinetic-Based Traffic Model 



Summarizing the above results, our gas-kinetic-based traffic model can be written in the general 
form (|l|), (0) of macroscopic traffic models, but with a non-local relaxation term. Since the 
equations are structurally related to, for example, the KKKL model, we find many similar non- 
linear phenomena p(J. This includes the sequence of stable, linearly unstable, and metastable 
regimes [|10| [58[] the local breakdown effect fl57| , the local cluster effect l lQfl , and, at sufficiently 
large densities, the formation of so-called dipole layers |58| . The response of equilibrium traffic 
to localized disturbances is similar to the Kerner-Konhauser model |57|. More specifically, for 
densities below a certain critical density p c i and above some density p C 4, homogeneous traffic is 
stable with respect to localized perturbations [Figure 0(a) + (d)], and for a range p C 2 < P < Pc3 of 
intermediate densities, it is linearly unstable, giving rise to cascades of traffic jams ["stop-and-go 
traffic", cf. Figure 0(c)]. For the two density regimes p c i < p < p c2 and p c3 < p < p c4 between 
the stable and the linearly unstable regions, traffic is metastable, i.e., it behaves unstable 



with respect to perturbations exceeding a certain critical amplitude Ap cr (p) [Figure 0(b)], but 
otherwise stable. For the self-organized density pj am inside traffic jams we find a typical value 
Pjam > Pc4 P^ - Furthermore, there exists a range p cv < p < p c3 with p cv > p c2 , where traffic is 
linearly unstable, but convectively stable, i.e., all perturbations grow, but they are eventually 
convected away in upstream direction 



Insert Figure [| about here. 

In addition, we obtain that, in the unstable traffic regime, the resulting flow-density relation 
differs from the equilibrium one (lying below the latter). We also find that the outflow Q out 
from traffic jams is independent of the initial conditions and the spatially averaged density 
[Figure |](a)]. Furthermore, the dissolution velocity C of traffic jams varies only a little with 
density [Figure f§(b)]. 

Insert Figure [| about here. 

The main difference with respect to other macroscopic traffic models is the non-local char- 
acter of the braking term. Nevertheless, our model can still be rewritten in the form of flux 
equations, with a non-local source term [51, |52|]. For this reason, we can apply various standard 
methods of numerical integration. It turns out that the non-local term has similar smoothing 
properties like a viscosity term. However, its effect is anisotropic due to the anticipation be- 
havior in forward direction reflecting that vehicles mainly react on density or velocity gradients 
ahead of them, not behind them (compare Fig. with Fig. [I]). Compared to a viscosity term, 
the non-locality does not change the hyperbolic character of the partial differential equations 
to a parabolic one. Therefore, it has favourable properties with respect to numerical stability 
and integration speed |37|. As a consequence, our model allows a robust real-time simulation 



of freeway stretches up to several thousand kilometers on a usual PC. 

Insert Figure [| about here. 



2.3 Model Calibration and Validation 

It turns out that our model can be easily calibrated to the static and dynamic properties 
of traffic flow data by a systematic procedure. In spatially homogeneous traffic, the average 
velocity relaxes to 



with 



V(p) 



VM = g | -i 



Vn 





1 



/ t[1 - p(p)]x(p)pA(p) l \ 
This also determines the equilibrium traffic flow per lane by 

Q e (p) = P^c(p) , 



A( Pr _ 



(19) 



(20) 



(21) 



which is depicted in Fig. 0(b). In the limit of small densities we find V c (p) ~ Vq, while in 
the limit of high densities (i.e. for (1 — p/p max ) < 1 or K C Vq), we have V e (p) ~ V(p), i.e. 
V e (p) ~ (1/p— 1/Pmax)/^- Apparently, the equilibrium flow-density relation is only affected by 
the model parameters Vo, T, and p max . The desired velocity Vq can be determined by fitting 
flow-density data at low densities by a straight line pVb, while the safe time headway T and 



the maximum density p max are obtained by fitting the data at high densities by a straight line 
(1 — p/pmax)/T, which crosses the abszissa at p max with a slope of —1/ (p max T). 

The model parameters r and 7 influence the stability behavior. Figure ^| shows that an 
increased value of r leads both to an increased range of instability, and to increased amplitudes 
(Pjam — Pout) of traffic jams. Further simulations showed that higher values of 7 tend to increase 
the stability of traffic. This is plausible, since 7 reflects the anticipation of future velocity 
changes. 7 and r also determine the shape and width of the downstream and upstream fronts 
connecting free and congested states (see Figs. [I] and |5|). Since r and 7 weakly influence the 
outflow Q out from traffic jams, the calibration of T, r, and 7 is repeated recursively until an 
optimal fit of the model to the empirical data is achieved. 

Insert Figure about here. 

Because of their intuitive meaning, the plausible range of values for the model parameters 
is rather restricted. V is given by the average free speed. p max must be consistent with 
the average length of vehicles plus a reasonable bumper-to-bumper distance of about 1.5 m in 
standing traffic. T should be compatible with the average time headway kept in homogeneous 
congested traffic. Reasonable (average) time headways are in the range T G [1.0 s, 2.5 s]. Initial 
accelerations a max = Vq/t are typically in the range a max G [1 m/s 2 , 4 m/s 2 ], corresponding 
to r G [10 s, 40 s]. (Notice that, for Vq = 158 km/h, r would have the meaning of the average 
acceleration time from to 100 km/h.) Therefore, a relaxation time r ~ 35 s is reasonable for 
freeway traffic, whereas r is smaller for city traffic. Finally, the minimum anticipation of traffic 
is to the car in front, implying 7 > 1. 

The above procedure of parameter calibration has been applied to single-vehicle data of 
the Dutch freeway A9, leading to the parameter set shown in Table |l|. It turns out that all 
optimized parameters have realistic values. In particular, this holds for r. The typical values 
for the safe time headway of T = 1.8 s is consistent with the rule "distance (in m) should not be 
less than half the velocity (in km/h)" suggested by German road authorities. For other data, 
however, one often finds that a somewhat smaller time headway yields a better fit. 

Insert Table [J about here. 

Since the model parameters are meaningful, it is simple to model changes of the traffic dynamics 
caused by external effects like environmental influences (cf. Figure |^). For example, a speed 
limit would be considered by decreasing Vo- Bad weather conditions leading to more defensive 
driving would be characterized by an increased time headway T and a lower value of Vq (plus 
a reduction of 7, if there is heavy fog). In rush-hour traffic, it is plausible to assume a higher 
percentage of experienced drivers than in holiday traffic, which would correspond to a higher 
7. Effects like a varying distribution of vehicle types can be modelled as well. For example, a 
higher proportion of lorries (long vehicles) would lead to a decrease of V and p max , but also to 
an increased value of r. 

Insert Figure [7| about here. 



2.4 Treatment of Boundaries in Open Systems 



Notice that the left-hand sides of the continuity equation ([[]) and the velocity equation (0) 
constitute a hyperbolic set of partial differential equations. Hence, the direction of informa- 
tion flow due to convective and dispersive processes is given by the characteristic lines pD 
Specifically, the velocities of information flow are given by 



2p ov \ 



1 dP\ dP 
— H 7— + ^— with 
4p 2 \dV) dp 



P(p, V) = P 9(p, V) 



(22) 



For all physically reasonable parameter sets they are always positive. Information flow against 
the direction of motion of the simulated vehicles is taken into account by the non-local inter- 
action term. The consequences are: (i) There is always a transport of information in both 
directions, upstream and downstream, (ii) The most suitable differencing scheme for the nu- 
merical integration is the upwind scheme pT| . 

By choosing the upwind scheme, the downstream boundary conditions are only relevant 
for the non-local interaction term. This is in accordance with the boundary conditions for 
microscopic traffic simulations, where the first car at the downstream boundary always needs 
a "phantom" leading car for its update process. In particular, the microscopic downstream 
boundary conditions determine only the acceleration of this first car, but do not determine 
the flow of exiting vehicles. At the same time, there always has to be a defined inflow of new 
vehicles at the upstream boundary. The same is true for the simulation of non-local macroscopic 
models with the upwind scheme. Therefore, it is no problem to formulate downstream boundary 
conditions which are always satisfied. 

At the upstream boundary, unphysical situations can emerge when a backwards moving 
jam approaches the boundary, but the assumed boundary (in)flow is higher than the flow 
inside the traffic jam. Then, the equations are no longer well-posed in a mathematical sense. 
Such situations may occur, whenever there is a slight offset in the arrival times of simulated 
and real jams at the boundary. This means the boundary conditions have to be implemented in 
a way that, on the one hand, copes with unexpected or offset density waves and, on the other 
hand, returns to the exact solution as soon as possible. This can be done by an algorithm that 
implements homogeneous von Neumann boundary conditions instead of Dirichlet boundary 
conditions for density and flow, if (i) there is congested traffic in the vicinity of the boundary, 
and (ii) the imposed inflow is larger than the traffic flow of the congested state |37j . Condition 

(i) means that, near the boundary, the local group velocity v g = d(pV e (p))/dp for instantaneous 
velocity relaxation is negative, or 

p(5,t)>p m (23) 

(with some small 5 > and p m being the density for maximum equilibrium flow). Condition 

(ii) can be expressed by requiring 

Q(0,t) > Q(6,t). (24) 

This algorithm has been tested as follows: First, we simulated a circular road of 20 kilometers 
length. The simulated scenario contains a local dipole-like perturbation in the initial conditions 
which develops to a backwards moving cluster [see Fig. |8](a)]. Then, we used the time-dependent 
results at the locations x = 5 km and x = 15 km as Dirichlet boundary conditions for the 
simulation of a new, 10 kilometer long open system. Figure |8](b) shows a comparison of the 
density profiles of the two systems at two times ti and t 2 . At ti = 8min, the cluster is not 
yet fully developed, and at t 2 = 80 min, it has propagated around the circle once. The profiles 
belonging to the same time points agree so well, that they would lay one upon each other. For 
this reason, we have depicted them with a slight artificial offset. After this hard and successful 
test, we are prepated to simulate the traffic phenomena observed on real freeways, which are 
open systems. 

Insert Figure || about here. 



2.5 Simulation of "Synchronized" Congested Traffic 



Recently, Kerner and Rehborn || have presented experimental data indicating a first-order 
transition to "synchronized" congested traffic. Traffic data from several freeways in Germany 
llL H, the Netherlands [§, |54|, [55|], and the United States fl8], |H indicate that synchronized 



traffic is the most common form of congested traffic. Synchronized traffic typically occurs at 
on-ramps when vehicles are added to already busy "freeways" and has the following properties: 

(i) The dynamics of the average velocities on all lanes is highly correlated ("synchronized"). 

(ii) The corresponding points in the flow-density plane are widely scattered, (iii) Synchronized 
traffic is characterized by a low average velocity, but, in contrast to traffic jams, the associated 
traffic flow is rather high, (iv) The transition to synchronized traffic is often caused by a 
localized and short perturbation of traffic flow that starts downstream of the on-ramp and 
propagates upstream with a velocity of about —10 km/h. (v) As soon as the perturbation passes 
the on-ramp, it triggers synchronized traffic which spreads upstream in the course of time, (vi) 
Downstream, synchronized traffic eventually relaxes to free traffic, (vii) Synchronized traffic 
often persists over several hours, (viii) The transition from synchronized traffic to free traffic 
occurs at a lower density and higher average velocity than the inverse transition (hysteresis 
effect). 

We believe that property (i), i.e., the synchronization requires a strong coupling among 
lanes. When the street becomes crowded, the drivers are trying to avoid obstructions by 
changing lanes. They fill developing gaps in the other lane until it becomes as slow as the own 
lane [Fig. §(a)], thereby balancing each upcoming disequilibrium among lanes. While restoring 



the balance is related with a considerable number of lane changes |53| , the lane changing rate 
becomes rather low as soon as the synchronization among lanes is reached. This idea can be 
indirectly supported by empirical traffic data which show that, even in congested traffic, there 
can be a significant difference in the velocities of cars and lorries [Fig. [|(c), (d)], which can be 
only attributed to lane changing and overtaking. 

Insert Figure [I] about here. 



Property (ii) seems to be the result of the diversity of vehicles ||55|| , as demonstrated in 
Figure f|(b). The other characteristic properties of "synchronized" traffic can be reproduced by 
the effective one-lane version of the non-local, gas-kinetic-based traffic model described above. 
This applies to homogeneous forms of "synchronized" traffic as well as to non-stationary forms 
(see oscillating congested traffic, OCT, discussed below). The essential point is that, along 
on-ramps (or off-ramps), there is a source term in the continuity equation ([I]). With 

v = 5=£ , (25) 
nL 

we assume that this source term is given by the actually observed inflow Q imp > from (or 
outflow Q rm p < to) the ramp, divided by the merging length L and by the number n of lanes. 



The inflow has an upper limit that depends on the downstream flow on the main road [p3[ . If 
one likes to model the tendency that some on-ramp vehicles enter as soon as possible regardless 
of the length of the ramp, one replaces L by an effective ramp length L* < L. Moreover, if the 
average velocity V rmp of entering or leaving vehicles on the ramp is different from the average 
velocity V on the main road, this gives an additional contribution 

-{V rmp -V) (26) 
P 

to the right-hand side of the velocity equation (fj), corresponding to a source term vV rrap in 
the related differential equation for the traffic flow Q = pV on the freeway |5^| . However, 
usually one assumes V rmp ~ V. 

Our simulations reproduced the observations (iii) to (viii) very realistically. For V rmp = V, 
the parameter values V = 128 km/h, p max = 160 vehicles/km, T = 1.6 s, r = 31 s, 7 = 1.0, and 
the observed boundary flows, [Figure 0(c)] , the simulated velocities and flows are, apart from 



fluctuations, in good agreement with the empirical curves presented by Kerner and Rehborn 
(compare Ref. [0 with Ref. ||). 



Insert Figure [77] about h 



ere. 



Our results suggest the following interpretation of the phase transition to "synchronized" 
congested traffic: Initially, the homogeneous flow Qmain upstream of an on-ramp is stable, while 
the higher downstream flow Qdown = Qmain + Qrmp/^ is metastable. Without any disturbance of 
the main or ramp flow, free traffic flow will continue. However, a (positive or even a negative) 
perturbation in the on-ramp or the main flow will eventually cause a breakdown of velocity 



close to the on-ramp. Speaking in the terminology of Daganzo et al. [62], the perturbation 
will eventually "activate the bottleneck" associated with the on-ramp. In the following, we will 
explain and quantify the mechanism of this transition. 



According to Figure ITUl the perturbation triggers a stop-and-go wave, which travels down- 



stream as long as it is small and upstream as it becomes larger, as is known from "localized 
clusters" H57| . Now, assume the downstream front of the cluster would pass the on-ramp. Then, 
since Q ma m [Figure [l^(c)] is lower than the characteristic outflow Q out from a jam, the cluster 
would eventually vanish. However, during its lifetime, the cluster would continue to emit the 
flow Qout, leading downstream of the ramp to a flow Qout + Qrmp/n > Qmax- As a consequence, 
as soon as the perturbation reaches the on-ramp, it induces congested traffic with a standing 
downstream front just at the end of the ramp. With an observed outflow Qout < Qout from 
synchronized traffic (with Q ou t ~ Qout), the average flow upstream is given by 

Qsync Qout Qrmp/^ • \^) 

Hence, the congested region upstream of the on-ramp is growing until the flow Qmam from the 
main road drops below Q syn c- Only then, the congested region starts to "melt", and it can 
easily take an additional hour or so until traffic returns to free flow, even if the maximum flow 
Q max has never been exceeded by (Q ma m + Qvm P /n)\ 

Now, consider the density p sync defined by Q syn c = Qc(Psync) in the congested part of the 
equilibrium flow-density relation Q c (p) (i.e. behind its maximum). If homogeneous traffic is 
(meta-)stable at p sync , the on-ramp induces "synchronized" congested traffic. The restriction 
Qrmp < Qout/2 |k§ implies Q sync > Qout(l - ^) and p sync < p jam , so that "synchronized" 
congested flow is significantly higher than the flow inside traffic jams. 



2.6 Phase Diagram of Traffic States at Bottlenecks 

Close to bottlenecks by on-ramps, lane closures, gradients, or other effects (including building 
sites along the road, bad street conditions, restricted visibility, curves, accidents in the opposite 
lanes, and slow vehicles), there can develop a variety of different traffic states. Free traffic 
(FT) and the above discussed type of "synchronized" congested traffic, which we will also 
call "homogeneous congested traffic" (HCT), are only two of them. Other congested traffic 
states that can emerge are "oscillating congested traffic" (OCT), "triggered stop-and-go waves" 
(TSG), "moving localized clusters" (MLC) and "pinned (standing) localized clusters" (PLC) 
(see Fig. [TTj) . All these different congested states are observable in empirical traffic data |T3], ITT], 
||, 0, [18], |6"2"| , although the "synchronized" congested states have sometimes been interpreted 
differently in the beginning. 



Insert Figure [TTJ about here. 



In our traffic simulations, we found out that most of the congested traffic states occur below 
the capacity limit Q max , and they are triggered by perturbations in traffic flow. This is only 
possible since the outflow Q out from traffic jams is an additional, self-organized capacity limit, 
which is smaller than Q max (due to the larger time headways occuring in accelerating traffic). 
The same fact brings about that most breakdowns of traffic flow could be avoided by control 
measures which manage to suppress these perturbations. 

Our simulations also indicate that the conditions, under which the different traffic states 



appear, are clearly defined |55], |22|]. This is illustrated by the phase diagram in Figure |TT|(c), 
displaying the traffic states as a function of the main flow Q ma in and the on-ramp flow per 
freeway lane, Q r mp/n (which is a measure of the flow reduction by a bottleneck). It explains 
why usually not all possible congested traffic states are observed at every bottleneck, but 
only a typical, site-specific subset of them. The phase diagram also allows to understand 
the various observed transitions between the traffic states, which correspond to crossing the 
phase boundaries separating the different states. Notice that these phase boundaries can be 
calculated analytically in terms of the flows Qd = Q e (p C i) characterizing the stability diagram 
(see Sec. |2.2|) . This implies that the general structure of the phase diagram in Fig. [ll|(c) is 
expected to be universal for all microscopic and macroscopic traffic models which possess the 
same instability diagram (i.e. the sequence of stable, metastable, unstable, convectively stable, 
metastable, and stable traffic with increasing density). Only the exact position and the shape 
of the phase boundaries will vary from one model to another. 



Multistability and Spatial Coexistence of States. Note that the phase diagram in 
Fig. 0(c) is for fully developed perturbations. For smaller perturbations, the phase boundaries 
will be shifted. The MLC and PLC regimes may even disappear, leading to the boundaries 
FT<->TSG, FT<->OCT, and FT^HCT. In any case, there are regions of multistability, where 

6§ 



we could 



the resulting traffic state depends on the initial condition. Similar to Lee et al. 
even find a small region of tristable traffic (TRI), where we can have FT, PLC, or OCT, depend- 
ing on the the respective "history" of the traffic dynamics. Figure |l2|(a) shows, for example, a 
transition from PLC to OCT. 

Another interesting point is the possible spatial coexistence of different traffic states. If 
congested traffic is linearly unstable, but convectively stable, upstream of a bottleneck we can 
have the following sequences of states, depending on the model parameters and the length of the 
congested region: (i) Empirical data and simulation results of the sequence HCT^OCT (with- 
out mergers of small density oscillations to fully developed stop-and-go waves) are presented 
in Ref. |7|. (ii) The sequence HCT^OCT^TSG (see Fig. 0(b)) is discussed in Ref. |l[ 
(for a slightly modified GKT model with frustration effects), and we think that it can explain 
the empirically observed sequence "synchronized" traffic^ "pinch region" — >stop-and-go waves 



reported in Ref. 



In Sec. |3.7|, we will come back to the subjects of this paragraph. 



Insert Figure 71 about here. 



Qualitative Interpretation of the Transitions. Finally, we give a qualitative interpreta- 
tion of the transitions between the various congested traffic states and the respective conditions, 
under which they appear. For example, the transition HCT-^OCT from homogeneous to oscil- 
lating congested traffic occurs, when the ramp flow Q rmp is decreased so that the synchronized 
flow (|27| ) falls into the unstable density regime instead of the (met a-) stable one. A further 
decrease in the ramp flow will reduce the average density in the congested region so much 
that we will have an alternation between regimes of free and congested traffic, which defines 
stop-and-go waves. We call them triggered stop-and-go waves (TSG) because, whenever an 



upstream travelling jam passes the bottleneck, it triggers a small perturbation which travels 
downstream as long as it is small, but eventually changes its speed as it grows and changes 
its propagation direction. Finally, the developed perturbation travels upstream and passes the 
bottleneck, where it gives birth to a new perturbation, and so on. 

Obviously, the small triggered perturbation cannot grow, when the traffic flow downstream 
of the bottleneck is (met a-) stable, which will lead to a single localized cluster (LC). The localized 
cluster will pass the bottleneck in upstream direction (MLC), if traffic flow is unstable or 
metastable, there. In contrast, the localized cluster will be pinned at the bottleneck (PLC), 
if traffic in the upstream region is stable, since a density cluster could not survive, there. If 
the traffic flow downstream of the bottleneck is also stable, only free traffic (FT) can persist. 
Finally, the pinned localized cluster will grow to an extended region of congested traffic, if its 
outflow Q out is smaller than the total flow (Qmain + Qrm P / n ) at the bottleneck. 



Discussion. Although there are also other proposals for an explanation of the various tran- 
sitions between free traffic, congested traffic, and stop-and-go waves |}I9|| , we see the following 
advantages of our approach: The transitions result naturally from a model for traffic flow on 
homogeneous roads, just by adding source terms (in case of on-ramp flows) or other kinds of 
inhomogeneities like speed limits (see Sec. |3.7|) , which are known to exist. The implementation 
of these inhomogeneities is straightforward, without the requirement of additional refined model 
ingredients. Moreover, analogous to the other parameters of the model, the inhomogeneities 
are easily measurable quantities. There is no model ingredient, which could not be relatively 
easily be verified or falsified by empirical studies. Furthermore, empirical results confirming the 



existence of a phase diagram are already available |67|. Finally, we think that the simulated 
traffic states related to inhomogeneities of the road arise so naturally, that any explanation of 
empirical data must take these states into account. 



2.7 Velocity Correlations of Successive Cars 

With a few exceptions only (see, e.g., Refs. [03, PJ), Boltzmann-like traffic equations are 



evaluated with the assumption of "vehicular chaos", which neglects correlations r between 
the velocities of successive vehicles. This is also true for the results presented above, mainly, 
because empirical data of the correlation coefficient r have not been available until very recently 
(see Fig. [13]). As expected, the velocity correlations are small for low vehicle densities, but 
considerable in congested traffic. Fortunately, it is possible to generalize the above non-local, 
gas-kinetic-based traffic model in a way that allows to take the velocity correlations into account 
[BBI, . It turns out that one has only to replace formula (Bf) by 



S = 9-2rV99 1 +9'. (28) 

Obviously, this mainly modifies the equilibrium velocity a bit, similar to a reduction of velocity 
variances. Therefore, we expect only quantitative, but no qualitative (fundamental) changes in 
the simulation results for finite values of r. 

Insert Figure [731 about here. 



3 Car-Following Models 

For about fifty years, now, researchers model freeway traffic by means of continuous-in-time 
microscopic models (car- following models) ||69|| . Since then, a multitude of car-following models 
have been proposed, both for single-lane and multi-lane traffic including lane changes. 



In the simplest case, the acceleration of an individual vehicle depends only on the distance 
to the vehicle in front. Well-known models of this type include the model by Newell [[70], or the 
"optimal- velocity model" by Bando et al. [7I|. To achieve a better anticipative driver behavior 
and to avoid collisions, the acceleration in other models depends also on the velocity and on 
the approaching rate to the front vehicle [73, |73|, [7JJ. Besides these simple models intended for 
basic investigations, there are also highly complex "high-fidelity models" like the Wiedemann 
model J75[ or MITSIM J7B| . They try to reproduce traffic as realistically as possible, but at the 
cost of plenty of parameters, which are difficult to calibrate. 

In the following, we will propose and discuss a novel follow-the-leader model that contains a 
few parameters only, all of which have a reasonable interpretation, are known to be relevant, are 
empirically measurable, and have the expected order of magnitude. The fundamental diagram 
and the stability properties of the model can be easily (and separately) calibrated to empirical 
data, and simulations of all kinds of observed traffic breakdowns are possible with the measured 
boundary conditions [B7|! Moreover, an equivalent macroscopic version of the model is known, 
which is not the case for most other microscopic traffic models. Finally, we mention that the 
new model behaves accident-free because of the dependence on the relative velocity. For similar 
reasons and because of metastability, it shows the self-organized characteristic traffic constants 
mentioned before. 



3.1 The Microscopic Intelligent-Driver Model (IDM) 



The acceleration assumed in the IDM is a continuous function of the velocity v a , the (netto) 
gap s a , and the velocity difference (approaching rate) Av a of vehicle a to the leading vehicle: 
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This expression is a superposition of the acceleration a^[l — {v a /v^) 5 } on a free road, and 
a braking deceleration Av a )/s a } 2 , describing the interactions with other vehicles. 

The deceleration term depends on the ratio between the "desired gap" s* and the actual gap 



s a , where the desired gap 
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is dynamically varying with the velocity and the approaching rate, reflecting an intelligent driver 
behavior. The IDM parameters are the desired velocity Vq, safe time headway T, maximum 
acceleration a, comfortable deceleration b, acceleration exponent 5, and the jam distances sq 
and s\. Furthermore, the vehicles have a finite length I which, however, has no dynamical 
influence. In Sections 



J] to |3.7] , we will assume identical "cars" , while in Sec. |3.8| we assume 
two different types, "cars" and "lorries" , see Table 0. For better readability, we will drop the 
vehicle index a in the following discussion of the model. 

Insert Table [| about here. 



3.2 Equilibrium Traffic of Identical Vehicles 

In equilibrium traffic (v a = 0, Av a = 0), drivers tend to keep a velocity-dependent equilibrium 
gap s e (v a ) to the front vehicle given by 
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(31) 



In particular, the equilibrium gap of homogeneous congested traffic (v <C vq) is essentially equal 

to the desired gap, s e (v) ~ s*(v,0) = so + sxJvJvq + vT, i.e., it is composed of a (small) 
high-density contribution, and a contribution vT corresponding to a time headway T. 

Solving Eq. ( |3ll) for the equilibrium velocity v = v e leads to simple expressions only for 
si = and 5 = 1, 5 = 2, or 5 — > oo. In particular, the equilibium velocity for the special case 
5=1 and so = s\ = is 

<^ = ^(- l + f^)- (32) 

Macroscopically, homogeneous traffic consisting of identical vehicles can be characterized by 
the equilibrium traffic flow Q e (p) = pV c (p) as a function of the traffic density p. For 5 = 1 and 
So = Si = 0, this "fundamental diagram" follows from Eq. ([32] ) together with the micro-macro 
relation between gap and density: 

1,11 
s = --l = . (33) 

P P Pmax 

The resulting relation for the density-dependent equilibrium velocity is identical with that of 
the GKT model, if the GKT parameter AA is set to zero [see Eqs. flT9| ) and (pOD I, which is a 
necessary condition for a micro-macro correspondence. Figure |l4](a) shows that the acceleration 
coefficient 5 influences the transition region between the free and congested regimes. For 5 — > oo 
and si = 0, the fundamental diagram becomes triangular-shaped: Q c (p) = min(f p, [1 — p(l + 
Sq)]/T). For decreasing 5, it becomes smoother and smoother. 

Insert Figure about here. 



3.3 Dynamic Single- Vehicle Properties 

Figures 0(b) and 0(c) show how the parameters a and b determine the acceleration and 
braking behavior of single vehicles. At t = 0, one vehicle starts with zero velocity and 2.5 
km of free road ahead. Initially, it accelerates with a and approaches smoothly the desired 
velocity vq. At x = 2.5 km, we assume a standing obstacle, for example, the end of a traffic 
jam. When approaching the obstacle, the model mimics "intelligent" drivers who anticipate 
necessary braking decelerations and brake so as not to exceed the comfortable deceleration b in 



normal situations [El 



3.4 Calibration 

The fundamental relations of homogeneous traffic are calibrated with vq (low density), 5 (tran- 
sition region), T (high density), and so and s\ (jammed traffic). The stability behavior of traffic 
in the IDM model is determined mainly by the model parameters a, b, and T. The density in 
and the outflow from traffic jams are also influenced by Sq and s±. Since the accelerations a and 
b do not influence the fundamental diagram, the model can be calibrated essentially indepen- 
dently with respect to traffic flows and stability. As in the GKT model, traffic becomes more 
unstable for decreasing a (which corresponds to an increased acceleration time r = v /a), and 
for decreasing T (corresponding to reduced safe time headways). Furthermore, the instability 
increases with growing b. This is also plausible, because an increased desired deceleration b 
corresponds to a less anticipative braking behavior. 



3.5 Collective Behavior and Stability Diagram 



Although we are interested in realistic open systems, it turned out that many features can be 
explained in terms of the stability behavior in a closed system. Figure ^(a) shows the stability 
diagram of homogeneous traffic on a circular road. The control parameter is the homogeneous 
(average) density p. We applied both a very small and a large localized perturbation to check 
for linear and nonlinear stability, and plotted the resulting minimum (p ou t) an d maximum 
(Pjam) densities after a stationary situation was reached. The resulting diagram is very similar 



to that of the macroscopic KKKL and GKT models |TD], 120]. In particular, it displays the 



following realistic features: (i) Traffic is stable for very low and high densities, but unstable for 
intermediate densities, (ii) There is a density range p c \ < p < p C 2 of metastability, i.e., only 
perturbations of sufficiently large amplitudes grow, while smaller perturbations disappear. Note 
that, for most IDM parameter sets, there is no second metastable range at higher densities, 
in contrast to the GKT and KKKL models, (iii) The density inside of traffic jams and the 
associated flow Qj am = Qc(Pjam), cf. Fig. |l5|(b), do not depend on p. As further "traffic 
constants", at least in the density range 20 veh./km < p < 40 veh./km, we observe a constant 
outflow (5 ou t = Qc(Pout) and propagation velocity v g = (Q ont - Qjam)/(Pout - Pjam) ~ -15 
km/h of jams. Figure |T5|(b) shows the stability diagram for the flows. In particular, we have 
Qci < Qout < Qc2, where Q ci = Q e (p C i), i>e., the outflow from congested traffic is metastable, 
while in the GKT, Q out m Q C 2 is only marginally stable. 

Insert Figure [73| about here. 

In open systems, a third type of stability becomes relevant. Traffic is convectively stable, 
if, after a sufficiently long time, all perturbations are convected out of the system. Both in 
the macroscopic models and in the IDM, there is a considerable density region p cv < p < p c3 , 
where traffic is linearly unstable but convectively stable. 

3.6 Microscopic Implementation of Bottlenecks 

In macroscopic simulations, a natural implementation of road inhomogeneities is given by on- 
and off- ramps, which appear as a source term in the continuity equation for the density. An 
explicit microscopic modelling of ramps, however, would require a multi-lane model with explicit 
simulation of lane changes. In order to avoid the associated complications, one can either apply 
the micro- macro link and simulate the ramp section macroscopically (see Sec. §), or introduce 
flow- conserving inhomogeneities by making one or more model parameters dependent on the 



location x of the road. Suitable parameters for the IDM are vo or T |2l], g7[. Local parameter 
variations act as a bottleneck, if the outflow Q' ont from congested traffic in the downstream 
section is reduced with respect to the outflow Q ont in the upstream section. This requires a 
reduced desired velocity v' < v or increased time headway T' > T, or both. Figure 0(a) 
shows a traffic breakdown induced by a linear decrease of the desired velocity from vq = 120 
km/h for x < —L/2 to v' = 95 km/h for x > L/2, while Fig. |H](c) shows the same effect for 
an increase of the safe time headway from T = 1.2 s to T' — 1.45 s (L = 200 m). Both flow- 
conserving bottlenecks result in a similar traffic dynamics which, however, depends strongly 
on the amplitude of the parameter variation. Qualitatively the same dynamics is observed in 
real traffic data [Fig. [16|(e)1 and in macroscopic models including on-ramps with a ramp flow 
of Qrmp/n = (Qout — Qout)- ^ ms suggests to define a general "bottleneck strength" AQ by 

AQ := Q mp /n + Q out - Q' out . (34) 

In particular, we have AQ = Q riap /n for on-ramp bottlenecks, and AQ = Q out — Q' ont for flow- 
conserving bottlenecks. In the following, we will vary Vq. The regions with locally decreased 



desired velocity can be interpreted as sections with uphill gradients (which reduce the maximum 
velocities of vehicles) . 

Insert Figure [T^ about here. 
3.7 Phase Diagram of Congested Traffic 

In contrast to closed systems, in which the long-term behavior and stability is essentially de- 
termined by the average traffic density, the dynamics of open systems is controlled by the 
inflow Qmain- Furthermore, traffic congestions depend on road inhomogeneities and, because 
of hysteresis effects, on the history of previous perturbations. For a given history, the traffic 
states can be summarized by a phase diagram spanned by Q ma in and AQ. Figure |17| shows 
the IDM phase diagram for traffic states that develop after a single density cluster crosses the 
inhomogeneity. Depending on Qmain and AQ, the initial perturbation (i) dissipates, result- 
ing in free traffic (FT), (ii) travels through the inhomogeneity as a moving localized cluster 
(MLC) and neither dissipates nor triggers new breakdowns, (iii) triggers a traffic breakdown to 
a pinned localized cluster (PLC), which remains localized near the inhomogeneity for all times 
and either is stationary, cf. Fig. |18|(a) for t < 0.2 h, or oscillatory (OPLC). (iv) Finally, the 
initial perturbation can induce extended congested traffic (CT), whose downstream boundaries 
are fixed at the inhomogeneity, while the upstream front propagates further upstream in the 
course of time. This kind of congested traffic can be homogeneous (HCT), oscillatory (OCT), or 
consist of triggered stop-and-go waves (TSG). The observed traffic states and their qualitative 
dependence on Q ma m and AQ is the same as for the GKT model, Fig. [II] when identifying AQ 
with Q T mp/ n according to Eq. fl3ip. In contrast to OCT, where there is permanently congested 
traffic at the inhomogeneity ("pinch region" [19, |2"]|), the TSG state is characterized by a series 



of isolated density clusters, each of which triggers a new cluster as it passes the inhomogeneity. 

Insert Figure [7?| about here. 

Boundaries between and Coexistence of Traffic States. Simulations show that the 
outflow Q ou t from the nearly stationary downstream fronts of OCT and HCT satisfies Q ont < 
Q' ont , where Q' ont is the outflow from clusters in homogeneous systems for the downstream model 
parameters. If the bottleneck is not too strong, we have Q oll t ~ Q'ont- Then, for all types of 
bottlenecks, the congested traffic flow is given by Q CO ng = Qout — Qrmp/ n ~ Q' ou t — Qrm P / n , or 

Qcong ~ Qont AQ , (35) 

which generalizes (|27|). Extended congested traffic (CT) only persists, if the inflow exceeds the 
congested traffic flow. Otherwise, it dissolves to PLC. This gives the boundary 

CT -> PLC : AQ « Q out - Q main . (36) 

If congested traffic flow is convectively unstable, the resulting oscillations lead to TSG or OCT. If 
it is linearly stable, Qcong < Qc3, we have HCT. If it is convectively stable, but linearly unstable, 
Qcong £ [Qc3,Qcv], one has a spatial coexistence of states with HCT near the bottleneck and 
OCT further upstream [21] (Fig. |18|), which is frequently found in empirical data of congested 



traffic. In the IDM, this frequent occurrence is reflected by the wide range of flows falling 
into this regime. For the "car" parameters, we have Q c ^ = 600 vehicles/h and Q cv = 1340 
vehicles/h. 



Insert Figure about here. 



Pinch Effect and Merging of Clusters. Careful investigations of traffic data related to 
OCT states showed two phenomena: (i) In a narrow region near the inhomogeneity, con- 



gested traffic is nearly stationary ("pinch region"), while further upstream, there are oscillations 
of the traffic density, (ii) While propagating upstream, the oscillations grow and merge to a few 
large-amplitude density clusters with free traffic in between. Detector data of other freeways, 
however, show a pinch effect without mergers, cf. Fig. |16](e). Simulating the IDM with the 
"car" parameters leads to very few mergers, cf. Figs. [TE|(c) and [T7]. For other parameters, 



however, the IDM reproduces mergers ending up with stop-and-go traffic ||21|| . A possible ex- 
planation is the "starvation effect": For the parameters chosen in this article, the outflow Q on t 
from density clusters is in the middle of the metastable region [Fig. |T3(b)], so medium-sized 



and large density clusters persist. For the parameters of Ref. ||21|| , however, the outflow from 
density clusters satisfies Q out ~ Q&, so only large-amplitude clusters survive, while all others 
dissipate. 

Multistability. In general, the local phase transitions between free traffic, pinned localized 
states, and extended congested states are hysteretic. In the regions between the two dotted lines 
of the phase diagram in Fig. [17|, both, free and congested traffic is possible, depending on the 
previous history. In particular, for all five indicated phase points (but not for the simulations 
of Fig. [16]), free traffic would persist without the downstream perturbation. In contrast, the 
transitions PLC — >OPLC, and HCT — >OCT — >TSG seem to be non-hysteretic, i.e., the type 
of pinned localized cluster or of extended congested traffic, is uniquely determined by Q ma i n 
and AQ. 



In a small subset of the metastable region labelled "TRI" in Fig. |17], we even found tristabil- 
ity between FT, PLC, and OCT. We obtained qualitatively the same also for the GKT model, 
and it has been found for the KKKL model with OPLC instead of PLC for the pinned localized 



state 0. Fi gure |l8](a) shows that a single moving localized cluster passing the inhomogeneity 



triggers a transition from PLC to OCT. Starting with free traffic, the same perturbation would 
trigger OCT as well, while we never found reverse transitions OCT — > PLC or OCT — > FT 
(without a reduction of the inflow). That is, FT and PLC are metastable in the tristable region, 
while OCT is stable. 



3.8 Multi-Species Single-Lane Traffic 

In this section, we assume heterogeneous single-lane traffic consisting of 70% "cars" and 30% 
"lorries" , where the latter are characterized by a lower desired velocity, lower accelerations, and 
a larger safe time headway compared to cars (Table 0). Again, we simulate a traffic breakdown 
to HCT at a flow- conserving inhomogeneity, where the desired velocity of cars is reduced from 
120 km/h to 65 km/h, and that of lorries from 80 km/h to 36 km/h. (Similar results are found 
for less dramatic increases of T.) To compare the result with real traffic data, we implement 
"virtual" detectors at several fixed locations. The detectors record passage times and velocities 
of each vehicle to determine the macroscopic flow Q = n T /r (where n T is the number of passing 
vehicles in the averaging interval r = 45 s), the arithmetic velocity average V, and the density 
P = Q/V- Figure |T9| shows the resulting fundamental diagram 3 km upstream of the bottleneck, 
and time series of Q and V at three upstream locations. For comparison, Fig. |l9|(d) and (e) show 
45 s averages of real single-vehicle data of the Dutch freeway A9 from Haarlem to Amsterdam 
on October 14, 1994. The detector is about 0.7 km upstream of the on-ramp causing the 
traffic breakdown. The simulated and real traffic data agree qualitatively, in particular in the 
following respects: (i) There is a wide scattering of flow-density data in the congested regime 
(looking like an anisotropic two-dimensional random walk), while the data occupy a nearly 



one-dimensional region in the free regime, (ii) The distribution of flow-density data shows the 
typical inverse-A form with a distinct gap between free and congested traffic data, (iii) During 
the breakdown, the velocity drops to 10-20 km/h, while the flow is reduced by only about 20%. 
(iv) In all regions, the relative fluctuations of the velocity are smaller than those of the flow. 
Near the bottleneck, the fluctuations of congested traffic flow are much smaller than those of 
free traffic, while further upstream, the fluctuations grow. 

Insert Figure \TQ about here. 

Notice that the fluctuations of the simulated data come essentially from the different ve- 
hicle types and not from deterministic instabilities. Macroscopically, the situation of Fig. |T^ 
corresponds to a spatial coexistence of HCT and OCT. (For a pure OCT or TSG state, there 
would be no gap between the flow-density data of free and congested traffic, both in real traffic 
data and in simulations f6?fl.) 



4 Micro- Macro Link 

One of the biggest problems in comparing microscopic and macroscopic models via coarse grain- 
ing is that a separation between the microscopic scale (of single vehicles) and the macroscopic 
scale (allowing to determine densities and average velocities) is hardly possible. Either the 
averaging intervals are too small and aggregate quantities like the density cannot be defined 
consistently, or the averaging intervals are too large so that the dynamics of the model is 
smoothed out. 



Generalizing an idea sketched in Ref. |5T|, we will now propose a way of obtaining macro- 
scopic from microscopic traffic models, which is different from the gas-kinetic approach and (at 
least) applicable to the case of identical driver-vehicle units. For this, we define an average 
velocity by linear interpolation between the velocities of the single vehicles a: 

y, t s = v a (t)[x Q -i(t) -x} + v a -i{t)[x - x a (t)] ^ 

Xa-xit) - X a (t) 

where x a _x > x > x a , and the index a of driver- vehicle units increases against the driving 
direction. While the derivative with respect to x gives us 

dV 

-7^ = [Va-l(t) - U a (t)]/[Xa-l(t) " X a (t)] , (38) 

the derivative with respect to t gives us the exact equation 

(s + 4) r =^"> (39) 

where 

_ a a {t)[x a -i(t) -x} + a a -i{t)[x - x a (t)] 

A[X ' ^ ~ *«-!(*)-**(*) (40) 

is the linear interpolation of the single-vehicle accelerations a a characterizing the microscopic 
model. For most car-following models, the acceleration function can be written in the form a a = 
Omic^a) Av a , s a ), where Av a = (v a — f a _i) is the approaching rate, and s a = (x a _i — x a — l a ) 
the (netto) distance to the vehicle in front. In order to obtain a macroscopic system of partial 
differential equations for the average velocity and density, the arguments v a , Av a and s a of 
the single-vehicle acceleration have to be expressed in terms of macroscopic fields as well. 



Specifically, we make the following approximations: A(x,t) ~ a m i C (V(x, t), AV(x, t), S(x, t)), 
with AV(x, t) = [V(x,t) — V a (x, t)], S(x, t) = |[p _1 (x, t) + p~ l (x,i)} — p max , and the non- 
locality given by g a (x, t) = g(x + l/p(x,t),t) with g G {p, V}. In this way we obtain a 
non-local macroscopic velocity equation, which supplements the continuity equation (P for the 
vehicle density and defines, for a given microscopic traffic model, a complementary macroscopic 
model. Note that it does not contain a pressure term (1/ p)dP '/ 'dx with P = p9, in contrast to 
the GKT model described above. 

As Fig. |2(] shows for the intelligent driver model, there are microscopic models for which 
the above procedure leads to a remarkable agreement. First, we have simulated an open free- 
way stretch with the macroscopic counterpart of the IDM with given macroscopic initial and 
boundary conditions. Then, we have simulated the first three kilometers of the stretch micro- 
scopically with the IDM [Fig. P0|(a), dark grey], using the same upstream boundary condition. 
As downstream boundary condition, we used the macroscopic simulation result at the interface 
(x = 3km) to the remaining macroscopic simulation section [Fig. |2D|(a), light grey]. The den- 
sity profile at the interface is illustrated in Figure |20|(b). Figure 0(c) shows the density profile 
1 km upstream of the interface obtained by the microscopic simulation, compared with that 
obtained from the separate macroscopic simulation of the whole system. The density p(x, t) in 
the microscopic simulation section for x a _i > x > x a was determined via the formula 

1 _ _L ^ g Q (t)[x a _i(t) -x] + s a -i(t)[x - x a (t)} 

p(x,t) (t)-X a (t) 
Insert Figure [2Z] about here. 

The micro-macro link of traffic models is only of practical relevance, if micro- and macrosim- 
ulations can be carried out simultaneously. For this, we have derived rules how to translate 
macroscopic initial and boundary conditions into microscopic ones, with almost identical sim- 
ulation results. In addition, we have developed methods to determine and apply the interface 
conditions (i.e. the respective upstream and downstream boundary conditions at interfaces) 
dynamically during the simulation, i.e. on-line. Since information must be able to flow through 
the interface in both directions, the formulation of dynamic interface conditions is a particu- 
larly tricky task which cannot be discussed, here. Figures ^l](a) and (b) show examples for 
the spatio-temporal evolution of the density on a circular road. One half is simulated with the 
microscopic IDM (dark grey) and the other half is simulated with the macroscopic counterpart 
of the same model (light grey). As can be seen, it is possible to connect both sections in a way 
that small perturbations propagating in forward direction as well as developed density clusters 
propagating backwards can pass the interfaces without any significant changes in the shape or 
propagation velocity. 



Insert Figure [||] about here. 



While traffic simulations with microscopic single-lane models are more intuitive and detailled 
than macroscopic simulations, a microscopic implementation of on- and off-ramps would require 
rather complex multi-lane models with lane-changing rules. Ramps can be treated much easier 
by macroscopic models, where they just enter by a simple source term in the continuity equation. 
Hence, the micro-macro link can be applied to combine both advantages. Figure shows an 
open system that is simulated essentially with a microscopic single-lane model (dark grey). 
Only the ramp section (light grey) is simulated macroscopically. The on-ramp flow triggers a 
traffic breakdown to a triggered stop-and-go state, i.e. a cascade of traffic jams which propagate 
in the upstream microscopic section. The qualitative features of this complicated dynamics are 
in very good agreement with a purely macroscopic simulation of the same system. 

Insert Figure 111 about here. 



5 Summary and Conclusions 



In this paper, we have shown that all presently known macroscopic phenomena of freeway 
traffic, including (i) the fundamental diagram, (ii) the characteristic parameters of congested 
traffic (outflow, propagation velocity, etc.), and (iii) the transitions between free traffic, "syn- 
chronized" or other congested traffic, and stop-and-go traffic can be reproduced and explained 
by microscopic and macroscopic traffic models based on plausible assumptions and realistic 
parameters, (iv) These models have also successfully reproduced a variety of empirical ob- 
servations in a semi-quantitative way, using the measured boundary conditions []55], 0. (v) 
Moreover, simply by simulating a mixture of different vehicle types, we obtained the observed 
scattering of congested flow-density data, although our traffic models are deterministic and have 



unique equilibrium relations, (vi) Generalizations to multi-lane traffic [53, 48] are expected to 



result in an even better agreement with empirical findings, e.g., a larger scattering of "virtual" 



detector data and flow-stabilizing effects of speed limits ]77|, |78| . (vii) Finally, we established 
a link between microscopic and macroscopic models, which allows to carry out simultaneous 
micro- and macrosimulations of neighboring freeway sections. To us, apart from developing 
measures for traffic optimization, the most interesting open question is which kind of observ- 
able phenomena are produced by the heterogeneity of driver-vehicle units in addition to the 
scattering of traffic data (see Ref. for an example). 

Acknowledgments 

The authors want to thank for financial support by the BMBF (research project SANDY, grant 
No. 13N7092) and by the DFG (grant no. He 2789). They are also grateful to Henk Taale and 
the Dutch Ministry of Transport, Public Works and Water Management for supplying the 
freeway data. 

References 

[1] W. Leutzbach, Introduction to the Theory of Traffic Flow, Springer, Berlin (1988). 
[2] D. Helbing, Verkehrsdynamik, Springer, Berlin (1997). 

[3] M. Schreckenberg and D. E. Wolf (eds.) Traffic and Granular Flow '97, Springer, Singapore 
(1998). 

[4] P. H. L. Bovy (ed.) Motorway Traffic Flow Analysis. New Methodologies and Recent Em- 
pirical Findings, Delft University Press, Delft (1998). 

[5] B. S. Kerner and H. Rehborn, Experimental properties of phase transitions in traffic flow, 
Physical Review Letters 49, 4030-4033 (1997). 

[6] D. Helbing and B. A. Huberman, Coherent moving states in highway traffic, Nature 396, 
738-740 (1998). 

[7] B. S. Kerner and P. Konhauser, Cluster effect in initially homogeneous traffic flow, Physical 
Review E 48, R2335-R2338 (1993). 

[8] R. D. Kuhne, Macroscopic freeway model for dense traffic — Stop-start waves and incident 
detection. In Proceedings of the 9th International Symposium on Transportation and Traffic 



Theory (edited by I. Volmuller and R. Hamerslag), pp. 21-42, VNU Science Press, Utrecht, 
The Netherlands (1984). 



[9] R. D. Kiihne, Freeway speed distribution and acceleration noise — Calculations from a 
stochastic continuum theory and comparison with measurements. In Proceedings of the 
10th International Symposium on Transportation and Traffic Theory (edited by N. H. 
Gartner and N. H. M. Wilson), pp. 119-137, Elsevier, New York (1987). 

B. S. Kerner and P. Konhauser, Structure and parameters of clusters in traffic flow, Physical 
Review E 50, 54-83 (1994). 

B. S. Kerner and H. Rehborn, Experimental properties of complexity in traffic flow, Phys- 
ical Review E 53, R4275-R4278 (1996). 

D. Helbing and M. Treiber, Gas-kinetic-based traffic model explaining observed hysteretic 
phase transition, Physical Review Letters 81, 3042-3045 (1998). 

B. S. Kerner and H. Rehborn, Experimental features and characteristics of traffic jams, 
Physical Review E 53, R1297-R1300 (1996). 

B. S. Kerner, Experimental characteristics of traffic flow for evaluation of traffic modelling. 
In Transportation Systems, Vol. II (edited by M. Papageorgiou and A. Pouliezos), pp. 793- 
798, International Federation of Automatic Control, Chania, Greece (1997). 

R. A. Raub and Pfefer, Vehicular flow past incidents involving lane blockage on urban 
roads: A preliminary exploration. In 77th TRB Annual Meeting, manuscript no. 004, 
Transportation Research Board, Washington, D.C. (1998). 

H. Y. Lee, H. W. Lee, and D. Kim, Empirical phase diagram of congested traffic flow, 
Preprint [http : / /xxx . lanl . gov/ abs/ cond-mat/9905292 . 

J. Treiterer and J. A. Myers, The hysteresis phenomenon in traffic flow. In Proc. 6th Int. 
Symp. on Transportation and Traffic Theory (edited by D. J. Buckley), pp. 13-38, Elsevier, 
New York (1974). 

C. F. Daganzo, M. J. Cassidy, and R. L. Bertini, Some traffic features at freeway bottle- 
necks, Transportation Research B 33, 25-42 (1999). 

B. S. Kerner, Experimental features of self-organization in traffic flow, Physical Review 
Letters 81, 3797-3800 (1998). 

M. Treiber, A. Hennecke, and D. Helbing, Derivation, properties, and simulation of a 
gas-kinetic-based, non-local traffic model, Physical Review E 59, 239-253 (1999). 

M. Treiber and D. Helbing, Explanation of observed features of self-organization in traffic 
flow. Preprint |http : //xxx . lanl . gov/abs/cond-mat/9901239| . 

D. Helbing, A. Hennecke, and M. Treiber, Phase diagram of traffic states in the presence 
of inhomogeneities. Physical Review Letters 82, 4360-4363 (1999). 

M. J. Lighthill and G. B. Whitham, On kinematic waves: II. A theory of traffic on long 
crowded roads, Proceedings of the Royal Society A 229, 317-345 (1955). 

P. I. Richards, Shock waves on the highway, Operations Research 4, 42-51 (1956). 



[25] G. F. Newell, A simplified theory of kinematic waves in highway traffic, Transportation 
Research B 27, 281-313 (1993). 

[26] C. F. Daganzo, The cell transmission model: A dynamic representation of highway traffic 
consistent with the hydrodynamic theory, Transportation Research B 28, 269-287 (1994). 

[27] D. F. Daganzo, The cell transmission model, Part II: Network traffic, Transportation Re- 
search B 29, 79-93 (1995). 

[28] M. Hilliges and W. Weidlich, A phenomenological model for dynamic traffic flow in net- 
works, Transportation Research B 29, 407-431 (1995). 

[29] J. P. Lebacque, A finite acceleration scheme for first order macroscopic traffic flow models. 
In Transportation Systems, Vol. II (edited by M. Papageorigiou and A. Pouliezos), pp. 
815-820, International Federation of Automatic Control, Chania, Greece (1997). 

[30] H. J. Payne, Models of freeway traffic and control. In Mathematical Models of Public 
Systems, Vol. 1 (edited by G. A. Bekey), pp. 51-61, Simulation Council, La Jolla, CA 
(1971). 

[31] M. Papageorgiou, Applications of Automatic Control Concepts to Traffic Flow Modeling 
and Control, Springer, Berlin (1983). 

[32] M. Cremer, Der Verkehrsflufi auf Schnellstrafien, Springer, Berlin, 1979. 

[33] W. F. Phillips, A kinetic model for traffic flow with continuum implications, Transportation 
Planning and Technology 5, 131-138 (1979). 

[34] H. Y. Lee, H. W. Lee, D., and Kim, Origin of synchronized traffic flow on highways and 
its dynamic phase transition, Phys. Rev. Lett. 81, 1130-1133 (1998). 

[35] G. B. Whitham, Linear and Nonlinear Waves, Wiley, New York (1974). 

[36] C. F. Daganzo, Requiem for second-order fluid approximations of traffic flow, Transporta- 
tion Research B 29, 277-286 (1995). 

[37] D. Helbing and M. Treiber, Numerical simulation of macroscopic traffic equations, Com- 
puting in Science and Engineering 1, 89-99 (1999). 

[38] M. Bando, K. Hasebe, K. Nakanishi, A. Nakayama, A. Shibata, and Y. Sugiyama, Phe- 
nomenological study of dynamical model of traffic flow, Journal de Physique I France 5, 
1389-1399 (1995). 

[39] D. Helbing and M. Schreckenberg, Cellular automata simulating experimental properties 
of traffic flow, Physical Review E 59, R2505-R2508 (1999). 

[40] D. Helbing, Gas-kinetic derivation of Navier-Stokes-like traffic equations, Physical Review 
E 53, 2366-2381 (1996). 

[41] D. Helbing, Derivation and empirical validation of a refined traffic flow model, Physica A 
233, 253-282 (1996). 

[42] I. Prigogine and R. Herman, Kinetic Theory of Vehicular Traffic, Elsevier, New York 
(1971). 



[43] S.-L. Paveri- Font ana, On Boltzmann-like treatments for traffic flow. A critical review of the 
basic model and an alternative proposal for dilute traffic analysis, Transportation Research 
9, 225-235 (1975). 

[44] P. Nelson, A kinetic model of vehicular traffic and its associated bimodal equilibrium 
solutions, Transport Theory and Statistical Physics 24, 383-409 (1995). 

[45] C. Wagner et al, Second order continuum traffic flow model, Physical Review E 54, 5073- 
5085 (1996). 

[46] A. Klar and R. Wegener, Enskog-like kinetic models for vehicular traffic, Journal of Sta- 
tistical Physics 87, 91-114 (1997). 

[47] P. Nelson and A. Sopasakis, The Prigogine- Herman kinetic model predicts widely scattered 
traffic flow data at high concentrations, Transportation Research B 32, 589-604 (1998). 

[48] S. Hoogendoorn and P. H. L. Bovy, A macroscopic multi-lane multi-class traffic flow model, 
Paper submitted to the 79th TRB Annual Meeting 1999 in Washington (1998). 

[49] A. Klar and R. Wegener, A hierarchy of models for multilane traffic I: Modeling, SIAM 
Journal on Applied Mathematics 59(3), 983-1001 (1999). 

[50] A. Klar and R. Wegener, A hierarchy of models for multilane traffic II: Numerical investi- 
gations, SIAM Journal on Applied Mathematics 59(3), 1002-1011 (1999). 

[51] D. Helbing, From microscopic to macroscopic traffic models. In A Perspective Look at 
Nonlinear Media. From Physics to Biology and Social Sciences (edited by J. Parisi, S. C. 
Miiller, and W. Zimmermann), pp. 122-139, Springer, Berlin, 1998. 

[52] D. Helbing and M. Treiber, Enskog equations for traffic flow evaluated up to Navier-Stokes 
order, Granular Matter 1, 21-31 (1998). 

[53] V. Shvetsov and D. Helbing, Macroscopic dynamics of multi-lane traffic, Physical Review 
E 59, 6328-6339 (1999). 

[54] D. Helbing, Empirical traffic data and their implications for traffic modeling, Physical 
Review E 55, R25-R28 (1997). 

[55] M. Treiber and D. Helbing, Macroscopic simulation of widely scattered synchronized traffic 
states, Journal of Physics A: Mathematical and General 32 L17-L23 (1999). 

[56] D. Helbing, Fundamentals of traffic flow, Physical Review E 55, 3735-3738 (1997). 

[57] B. S. Kerner, P. Konhauser, and M. Schilke, Deterministic spontaneous appearance of 
traffic jams in slightly inhomogeneous traffic flow, Physical Review E 51, R6243-R6246 
(1995). 

[58] B. S. Kerner, P. Konhauser, and M. Schilke, "Dipole- layer" effect in dense traffic flow, 
Physics Letters A 215, 45-56 (1996). 

[59] M. C. Cross and P. C. Hohenberg, Pattern formation outside of equilibrium, Review of 
Modern Physics 65, 851-1112 (1993). 

[60] W. F. Ames, Nonlinear Partial Differential Equations in Engineering, Vols. I & II, Aca- 
demic Press, New York (1965). 



[61] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in 
C: The Art of Scientific Computing, Cambridge University Press, Cambridge (1992). 

[62] C. F. Daganzo, M. J. Cassidy, and R. L. Bertini, Possible explanations of phase transitions 
in highway traffic, Transportation Research A 33, 365-379 (1999). 

[63] C. F. Daganzo, The nature of freeway gridlock and how to prevent it. In Proc. 13th 
Int. Symp. on Transportation and Traffic Theory (edited by J. B. Lesort), pp. 629-646, 
Pergamon-Elsevier, New York (1996). 

[64] D. Helbing, Modeling multi-lane traffic flow with queuing effects, Physica A 242, 175-194 
(1997). 

[65] D. Helbing and M. Treiber, Jams, waves, and clusters, Science 282, 2001-2003 (1998). 

[66] H. Y. Lee, H. W. Lee, and D. Kim, Dynamic states of a continuum traffic equation with 
on-ramp, Physical Review E 59, 5101-5111 (1999). 

[67] M. Treiber, A. Hennecke, and D. Helbing, Congested traffic states in empirical observa- 



tions and microscopic simulations, preprint |http : // arXiv . org/ abs/ cond-mat/0002177| , 
submitted to Physical Review E (2000). 

[68] D. Helbing, A. Hennecke, V. Shvetsov, and M. Treiber Macroscopic traffic simulation based 
on a gas-kinetic, non-local traffic model, Transportation Research B, in print (2000). 

[69] A. Reuschel, Fahrzeugbewegungen in der Kolonne, Osterreichisches Ingenieur-Archiv 4, 
193-215 (1950). 

[70] G. F. Newell, Nonlinear effects in the dynamics of car following, Operations Research 9, 
209-229 (1961). 

[71] M. Bando, K. Hasebe, A. Nakayama, A. Shibata, and Y. Sugiyama, Dynamical model of 
traffic congestion and numerical simulation, Physical Review E 51, 1035-1042 (1995). 

[72] P. G. Gipps, A behavioural car-following model for computer simulation, Transportation 
Research B 15, 105-111 (1981). 

[73] S. Kraufi, Microscopic Modelling of Traffic Flow, PhD thesis, DLR, FB 98-08, Cologne, 
(1998). 

[74] D. Helbing and B. Tilch, Generalized force model of traffic dynamics, Phys. Rev. E 58, 
133-138 (1998). 

[75] R. Wiedemann, Simulation des Strafienverkehrsflusses, Institut fur Verkehrswesen, Uni- 
versitat Karlsruhe (1974). 



[76] See internet page [http : / /hippo .mit . edu/products/mitsim/main . html 



[77] Interactive simulations of the multi-lane IDM are available at 
|http : //www.uni-stuttgart . de/treiber/MicroApplet/ . 



[78] R. Sollacher and H. Lenz, Nonlinear control of stop-and-go traffic. In Traffic and Granular 
Flow '99 (edited by D. Helbing, H. J. Herrmann, M. Schreckenberg, and D. E. Wolf), 
Springer, Berlin, 2000. 



Table 1: Typical values of model parameters used in the non-local, gas- kinetic-based traffic 
model to reproduce empirical data of the Dutch freeway A9. 
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Table 2: Model parameters of the IDM model used throughout this paper. 
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80 km/h 


1.7 s 
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10 m 
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Figure 1: Simulation of a downstream front with a jam density of 140 vehicles/km/lane. 
Shown is the spatio-temporal development (a) of the density p(x,t), and (b) the flow Q(x,t) = 
p(x,t)V(x,t). Free boundary conditions were used on both sides. 




Figure 2: Comparison of (a) the density-dependent relative variance in units of the squared 
average velocity, and (b) the equilibrium flow-density relation ([H]) of our gas-kinetic-based 
traffic model (solid lines) with empirical data (crosses). The empirical data were obtained from 
single- vehicle data of the Dutch motorway A9 (where a speed limit of 120 km/h applies), by 
averaging over one-minute intervals. Note that both, the variance-density relation and the 
velocity-density relation in equilibrium, are fitted by a single, density-dependent function A(p). 
It turns out that the deviation of the variance-density relation from the empirical data at small 
densities is not of great importance for the dynamics of the traffic model. It can, however, be 
accounted for by the function W(p), see formulas (|i"3|) and (|i"4|). 





Figure 3: Spatio-temporal evolution of the traffic density p(x,t) on a uni-directional ring of 
circumference 10 km, starting with homogeneous traffic to which a localized initial perturbation 
of amplitude 10 vehicles/km/lane is added, (a) Free and stable traffic at an average density 
of 15 vehicles/km/lane, (b) metastable traffic for 25 vehicles/km/lane, (c) unstable traffic for 
35 vehicles/km/lane, (d) stable congested traffic for 55 vehicles/km/lane. The chosen model 
parameters are displayed in Table |l|. 
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Figure 4: (a) Characteristic flows resulting in fully developed stop-and-go traffic on a circular 
road, as a function of the spatially averaged density. Depicted are the flows Qjam hi the 

jammed regions ( ), the outflows Q out from jams ( ), and the average flows ( — ). For 

comparison, the equilibrium flow Q e {p) = pV e (p) with V e (p) from equation ( |i9|) is also shown 
(•••). Notice that, in the unstable range, the average dynamic flow is lower than the equilibrium 
flow, (b) Propagation velocity of density clusters as a function of average density. We started 
our simulations with an initial perturbation of amplitude 10 vehicles/km/lane and measured 
the group velocity after a sufficiently long transient time. The propagation velocity at low 
densities is positive, but slower than the average vehicle velocity. In the instability region, the 
negative propagation velocity of fully developed traffic jams (i.e. their dissolution velocity) 
is independent of the initial density, and its magnitude is in agreement with empirical data. 
Notice that the dissolution velocity C of jam fronts corresponds to the slope of the average flow 
depicted in (a). 




Figure 5: Simulation of an upstream front with initial densities of p\ = 15 vehicles/km/lane 
and p2 = 140 vehicles/km/lane. Shown is the evolution of (a) the density p(x,t), and (b) the 
flow Q(x,t) = p(x,t)V(x,t). Note that, in contrast to other macroscopic traffic models fl36| , 
the simulations do not produce densities above p max or negative flows. In (b), the direction of 
the space and time axes is reversed for illustrative reasons. 
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Figure 6: Stability diagram for a localized perturbation of homogeneous traffic on a ring of 
circumference 10 km. Both diagrams show the developing maximum and minimum densities 
Pj am and p out as a function of the spatially averaged density, measured after a dynamical equi- 
librium was reached. The unstable traffic regime corresponds to the density range where the 
jam amplitude (pj am — Pout) is large (rectangle-like shaped regions). Diagram (a) shows the 
dependence of the stability diagram on the perturbation amplitude Ap. One can clearly see 
two density ranges [p c i,p C 2] with p c i = 21 vehicles/km/lane and p C 2 = 24 vehicles/km/lane, 
and [p C 3, Pd] with p c3 = 51 vehicles/km/lane and p C 4 = 55 vehicles/km/lane, where traffic is 
non-linearly stable, i.e., stable for small perturbations, but unstable for large perturbations. 
In the range [p C 2,Pc3], homogeneous traffic is unstable for arbitrary perturbation amplitudes. 
Diagram (b) shows the stability diagram for various relaxation times r and a perturbation 
amplitude of Ap = 1 vehicle/km/lane. 




Figure 7: Equilibrium flow-density relations of the non-local, gas-kinetic-based traffic model. 
The diagrams (a) and (b) show the variation with the model parameters Vq and T, respectively. 
Notice that a reasonable speed limit almost does not affect the flow in the congested regime, 
while a variation of the following time t can account for large variations in the flow in the 
congested density regime. In each diagram, the solid lines correspond to the standard parameter 
set displayed in Table p]. 




Figure 8: (a) Circular reference system simulated in order to obtain boundary conditions for the 
simulation of a comparable open system, (b) Comparison of density profiles at two times t\ = 
8 min and ti = 80 min, obtained in the cirular reference system (solid) and in an open system 
which was simulated with boundary conditions taken from the reference scenario (dashed). In 
order to make the curves distinguishable, one of the profiles was shifted by 100 m. 
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Figure 9: (a) Difference between the average velocity V^(p) in the left lane and Vi(p) in the 
right lane as a function of the lane-averaged vehicle density p. Points correspond to 1-minute 
data from November 2, 1994, which were determined from single-vehicle data of the Dutch 
freeway A9. We find no gap in the data cloud, so that there is no support for the hypothesis 
of a phase transition to synchronization among lanes. The same conclusion can be drawn from 
the averages of the 1-minute data belonging to the same density, which are represented by 
diamonds. We find a linear decrease without discontinuities up to a density of 40 vehicles per 
kilometer and lane, above which the velocity in both lanes is the same. 

(b) The displayed points in flow-density space correspond to lane-averaged one-minute data 
(dark crosses) and related simulation results (grey boxes) at the location of an on-ramp (af- 
ter |55|). The simulations manage to reproduce both, the quasi-linear flow-density relation at 
small densities and the scattering of the data over a two-dimensional region in so-called "syn- 
chronized" congested traffic at medium and high densities. The simulations are based on the 
assumption that, in an effective model of mixed traffic, the time-dependent parameters X(t) 
are given by a linear interpolation X(t) = {p\ (t)Xi + [1 — p\ (t)]X cai } between the parameters 
X car of cars and X\ Q of lorries, where p\ (t) denotes the time-dependent fraction of long vehicles 
passing the cross section. By lines, we have displayed the assumed equilibrium flow-density 

relations for traffic consisting of 100% cars ( — ), and of 100% long vehicles ( ). (c), (d) 

Average velocities of cars and lorries (long vehicles) in the left and the right lane as a function 
of density. The difference of these empirical curves shows a minimum around 25 vehicles per 
kilometer, where cars are almost as slow as lorries (after JJ). However, at higher densities, cars 
are faster again, which shows that there must be overtaking maneuvers at medium and high 
densities, which does not support the hypothesis of synchronization among lanes. 
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Figure 10: Spatio-temporal evolution of the lane-averaged density after a small peak of inflow 
from the on-ramp (after ||12|| ). The on-ramp merges with the main road at x — km with 
a merging length of L = 300m. Traffic flows from left to right. In (a), the parabolically 
shaped region of high density corresponds to "synchronized" congested traffic. Plot (b) shows 
the formation of this state in more detail. The time- dependent inflows Q ma i n at the upstream 
boundary and Q rmp /n at the on-ramp are displayed in (c). 




Figure 11: (a), (b), (d), (e) Spatio-temporal dynamics of typical representatives of congested 
states which are triggered at an on-ramp by a fully developed perturbation travelling upstream. 
The middle of the on-ramp is located at x — 8.0 km. The respective states are (a) triggered stop- 
and-go traffic (TSG) for Qmam = 1660 vehicles/h/lane and Q r m P /n = 75 vehicles/h, (b) oscil- 
latory congested traffic (OCT) for Q ma in = 1540 vehicles/h/lane and Q rmp /n = 170 vehicles/h, 
(d) a pinned localized cluster (PLC) for Qmain = 1450 vehicles/h/lane and Q rmp /n = 60 ve- 
hicles/h, and (e) homogeneous congested traffic (HCT) for Q mSLin = 1350 vehicles/h/lane and 
Qrmp/n = 400 vehicles/h. The other model parameters used in the above simulations are 
V = 110 km/h, p max = 140 Fz/km, r = 40 s, T = 1.7 s, 7 = 1.2, A = 0.008, AA = 0.02, 
p c = 0.27p max , Ap = 0.1p ma x, and L = 400 m. They can be considered as typical for German 
freeways, where all the above congested traffic states have been empirically observed, 
(c) Corresponding phase diagram of the traffic states forming in the vicinity of an on-ramp as 
a function of the inflows Qmain and Q rmp /n per freeway lane on the main road and the on-ramp 
for the above parameters. The different states are classified at t — 90 min, i.e., after a suffi- 
ciently long transient period. Displayed are homogeneous congested traffic (HCT), oscillatory 
congested traffic (OCT), triggered stop-and-go traffic (TSG), moving localized clusters (MLC), 
pinned localized clusters (PLC), and free traffic (FT). Spatially coexisting states are indicated 
with a plus sign. The broken line represents the condition Qmain + Qrm P / n = Qm&x that sepa- 
rates the region, in which a breakdown of traffic flow is caused by exceeding the theoretically 
possible freeway capacity Qmax (upper right corner). Note that all congested traffic states below 
this line are caused by perturbations and, therefore, avoidable by technical control measures. 




Figure 12: (a) Metastable pinned localized cluster (PLC) and stable oscillating congested 
traffic (OCT), simulated with the non-local GKT model (with model parameters vq = 120 
km/h, T = 1.8 s, r = 50 s, p max = 130 vehicles/km, 7 = 1.2, A = 0.008, and AA = 0.008). 
The inflow to the main road is 1400 vehicles per hour and lane, and the assumed inhomogeneity 
of traffic flow comes from an on-ramp of length 500 m with a constant inflow of 140 vehicles 
per hour and lane. The metastable PLC is triggered by a triangular-shaped density peak in 
the initial conditions (of total width 500 m, centered at x = 0), in which the density rises from 
14 vehicles/km to 40 vehicles/km. The OCT is triggered by a density wave introduced by the 
downstream boundary conditions, (b) Spatiotemporal evolution of the traffic density according 
to the gas-kinetic-based traffic model ||20|| . The assumed inhomogeneity of traffic flow comes 
from an on-ramp of length 200 m with an inflow of 220 vehicles per hour and freeway lane. 
The inflow to the main road is 1570 vehicles per hour and lane, and the breakdown of traffic 
flow is triggered by a perturbation AQ(t) of the inflow with a flow peak of 125 vehicles per 
hour and lane (see main text). The assumed model parameters are Vq = 120 km/h, T = 1.5 s, 
t = 30 s, p m ax = 120 vehicles/km, and 7 = 1.2, while the parameters for the variance prefactor 
[0 A(p) = A Q + AA{tanh[(p - p c )/Ap] + 1} are A = 0.008, AA = 0.02, p c = 0.27p max , 
and Ap = 0.1p max . In order to have a large region of linearly unstable but convectively stable 
traffic, we introduced a "resignation effect", i.e., a density-dependent reduction of the desired 
velocity V to V£{p) = V - AV/{1 + exp[(p^ - p)/Ap']} with AV = 0.9V , p' c = 0.45p max , and 
Ap = 0.1p max . 




Figure 13: Correlation between the velocity of successive vehicles (a) in the right lane as a 
function of the density (b) as a function of the gap between successive vehicles for a density 
of 10 vehicles per kilometer. The relations were determined from single-vehicle data from the 
Dutch freeway A9. According to our empirical results, the correlation coefficient is about zero 
at small densities and reaches a maximum at around 20 to 30 vehicles per kilometer, where the 
transition from free to congested traffic occurs. Afterwards, it stays on a high level (around 
0.65). In the left lane, the velocity correlation is a little bit higher, probably because of the 
smaller fraction of lorries. The correlation coefficient is not sensitive to the type of the leading 
vehicle (car or lorry), but to that of the following vehicle, which is reasonable. However, even 
at small densities, the correlation depends strongly on the gap to the next vehicle ahead and 
becomes almost one for very small gaps. This is, because different velocities would imply a 
high danger of accidents, then. 




Figure 14: (a) Equilibrium flow-density diagram of identical IDM vehicles with variable accel- 
eration exponent 5 and the "car" parameters of Table ^otherwise, (b), (c) Temporal evolution 
of the velocity and acceleration of a single vehicle approaching a standing obstacle, which is 
reached after t = 106 s (cf. main text). The IDM parameters are a = 1 m/s 2 , 6 = 2 m/s 2 , and 
the "car" parameters of Table ^| otherwise. 
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Figure 15: Stability diagram of homogeneous traffic (for "car" parameters) in a closed system as 
a function of the homogeneous density p for small (grey) and large (black) initial perturbations 
of the density In plot (a), the upper two lines display the density inside of density clusters 
after a stationary state has been reached. The lower two lines represent the density between 
the clusters. Plot (b) shows the corresponding flows and the equilibrium flow-density relation 
(dotted). The critical densities p C j and flows Qd are discussed in the main text. 




Figure 16: Traffic breakdowns at differently implemented bottlenecks, (a) IDM simulation of 
the spatio-temporal density (inflow Q ma i n = 1670 vehicles/h) for a bottleneck corresponding to a 
decrease of Vq in the downstream region. The density is derived from the microscopic distance s 
via relation ( |5BD and a subsequent linear interpolation, (b) Related equilibrium flows upstream 
and downstream, (c), (d) Bottleneck corresponding to an increase of the safe time headway T. 
(e) Density obtained from 1-min detector data of the German freeway A9-South on Oct. 29, 
1998. The traffic breakdown takes place upstream of the intersection "Neufahrn" at x — 512 
km. 




Figure 17: Phase diagram of congested traffic states and corresponding spatio-temporal density 
plots as simulated with the IDM. The control parameters are the inflow Q ma i n and the bottleneck 
strength AQ, which increases with (v - v' Q ). The traffic states FT, HCT, OCT, TSG, MLC, 
and (O)PLC are explained in the main text. "TRI" indicates a tristable region. 




Figure 18: (a) Transition from a pinned localized cluster (PLC) to oscillating congested traffic 
(OCT) in the tristable region (Qmain — 1480 vehicles/h, v' = 24 m/s corresponding to AQ = 
220 vehicles/h), triggered by a large perturbation, (b) Spatial coexistence of HCT and OCT 
for Q ma i n = 1350 vehicles/h and v' = 16 m/s corresponding to AQ = 400 vehicles/h. 
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Figure 19: (a) Flow-density diagram and (b),(c) time series of single-lane heterogeneous traffic. 
(d),(e) Empirical data of extended congested traffic on the Dutch freeway A9. 
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Figure 20: Results of simulations with the microscopic IDM (dark grey) and its macroscopic 
counterpart (light grey), assuming identical boundary and interface conditions. 




Figure 21: Results of simulations of a circular road, half of which was simulated microscopically 
(dark grey), while the other half was simulated macroscopically (light grey). Whereas the left 
picture shows the forward propagation of a decaying initial perturbation at small vehicle density, 
the right one illustrates the development of a backwards moving traffic jam at medium density. 




Figure 22: Result of a simultaneous micro-macro-simulation of triggered stop-and-go traffic, 
which is caused by high inflow from the ramp at x — 11km. While the ramp section (light 
grey) is implemented macroscopically, the other two sections (dark grey) were simulated mi- 
croscopically. 



